LMXBs may be important LIGO sources after all 
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ABSTRACT 

Andersson et al. and Bildsten proposed that the spin of accreting neutron 
stars is limited by removal of angular momentum by gravitational radiation which 
increases dramatically with the spin frequency of the star. Both Bildsten and 
Andersson et al. argued that the r— modes of the neutron star for sufficiently 
quickly rotating and hot neutron stars will grow due to the emission of gravita- 
tional radiation, thereby accounting for a time varying quadrupole component 
to the neutron star's mass distribution. However, Levin later argued that the 
equilibrium between spin-up due to accretion and spin-down due to gravitational 
radiation is unstable, because the growth rate of the r— modes and consequently 
the rate of gravitational wave emission is an increasing function of the core tem- 
perature of the star. The system executes a limit cycle, spinning up for several 
million years and spinning down in less than a year. However, the duration of 
the spin-down portion of the limit cycle depends sensitively on the amplitude at 
which the nonlinear coupling between different r— modes becomes important. As 
the duration of the spin-down portion increases the fraction of accreting neutron 
stars which may be emitting gravitational radiation increases while the peak flux 
in gravitational radiation decreases. Depending on the distribution of quickly 
rotating neutron stars in the Galaxy and beyond, the number of gravitational 
emitters detectable with LIGO may be large. 

Subject headings: stars: neutron — gravitational waves — stars: oscillations 

1. Introduction 

Accretion onto the surface of a neutron star can in principle spin up the rotation of the 
neutron star until the spin frequency equals the Kepler frequency of the inner edge of the disk. 
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In low-mass x-ray binaries, the disk is thought to extend to stellar surface so the maximal 
frequency that the neutron star can achieve exceeds 1 kHz. However, the observed and 
inferred spin frequencies of neutron stars in low-mass x-ray binaries (LMXBs) are clustered 
around 250-500 Hz (e.g Bildsten 1998); the millisecond X-ray pulsars SAX J1808.4 and 
XTE J1751 have slightly higher frequencies of 402 Hz and 435 Hz (Wijnands & van der 
Klis 1998; Markwardt & Swank 2002). Millisecond radio pulsars have been discovered with 
frequencies up to 640 Hz (Backer et al. 1982). All of these limits are well below the Keplerian 
limit on the spin frequency of a neutron star, so an alternative explanation for the maximal 
observed spin frequency of neutron stars is required. Andersson et al. (1999) and Bildsten 
(1998) proposed that inertial modes (specifically the r— modes) inside the neutron may grow 
while generating gravitational waves (GW). For sufficiently quickly rotating stars, GW carry 
away the angular momentum as quickly as it is deposited on the star by accretion. 

Levin (1999) found that this proposed equilibrium between spin-up and spin-down is 
unstable. A neutron star will execute a limit cycle (see Andersson et al. 2000, for additional 
discussion). It spins up for several million years and then quickly spins down emitting GW 
in less than a year. Only a small fraction of neutron stars is spinning down at any time; 
it is unlikely that any neutron star within the galaxy is currently spinning down, so none 
would be detected by LIGO. However, the duration of the spin down depends sensitively on 
the assumed maximal amplitude of the r— mode. Levin (1999) assume that the r— modes 
saturate when their amplitude is of order unity. Sperhake et al. (2001), Wu et al. (2001) 
and Arras et al. (2002) found that the saturation amplitude may be two to three orders of 
magnitude smaller; this may increase the duration of spin down to be greater than several 
thousand years. This dramatically increases the number of neutron stars whose GW could 
be detected both by increasing the typical GW flux relative to the estimates of Andersson 
et al. (1999) and Bildsten (1998) and the number of currently emitting sources relative to 
Levin (1999). 

This Letter explore the implications of saturation of r— modes in rapidly rotating neu- 
tron stars whose spins are accelerated by accretion. §2 will describe a series of straightforward 
calculations similar to those of Levin (1999) but with various values of the saturation ampli- 
tude (§2.1). Both Andersson et al. (2000) and Andersson & Kokkotas (2001) have estimated 
the number of observable sources if the duty cycle corresponds to an r— mode saturation 
amplitude of unity. The observed number of millisecond pulsars, the presumed descendents 
of LMXBS (e.g. Bhattacharya & van den Heuvel 1991), yields an estimate of the number of 
potential sources and their amplitudes (§§2.2-2.3) as a function of the duty cycle. §3 will 
outline some consequences of these results and the observational outlook. 
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2. Calculations 

2.1. Thermal and spin evolution 

The calculations presented here are essentially a recapitulation of those of Levin (1999); 
however, a brief summary of the equations governing the system is useful (Owen et al. 
1998). The angular velocity of the star il is normalized by the dynamical frequency of 
the star, (l = VL/ (nGp) 1 ^ 2 where p is the mean density of the star. If the dimensionless 
r— mode magnitude is less than the amplitude at which the mode saturates a max , the spin 
and r— mode amplitude are determined by 

dn 2a 2 Q Cl [AIM 

~dt ~ ~l + a 2 Q7 v + V3/M P ^' 
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dt \ r grav T v l + Q 2 Q / 

where Q and / are determined by the density profile of the star, a n = 1 polytrope (Lindblom 
et al. 1998). Following Levin (1999), Q = 0.094 and / = 0.261. p depends on the angular 
velocity at the inner edge of the accretion disk relative to that of the star. The time for 
the star to spin up is inversely proportional to p, so it can be subsumed into the uncertain 
value of the accretion rate M, i.e. p — 1. M is assumed to be 10 _8 M Q yr" 1 and M is 
the mass of the star. The timescale for the growth and attenuation of the r— mode for a 
fluid neutron-star consisting of neutrons, protons and electrons are r grav = — 3.26f2~ 6 s and 
T' 1 = (1.03 x lO^)- 1 ^ -2 + (6.99 x lO^s)- 1 ^ 6 where T 8 = T/(10 8 K) (Owen et al. 1998). 

Both the growth and decay rates depend on the properties of the neutron star. The 
presence of a crust (Lindblom et al. 2000; Wu et al. 2001; Levin & Ushomirsky 2001), 
magnetic field (Rezzolla et al. 2000), hyperonic or superfluid core (Lindblom & Mendell 
2000; Lindblom & Owen 2002) can dramatically affect the dynamics of r— modes in accreting 
neutron stars, shrinking the range of spin frequencies and temperatures where an r— mode 
will grow. On the other hand, a strange-quark-matter core may decrease the viscosity, 
shifting the instability region to cooler temperatures (Andersson et al. 2001). Consequently, 
the observation of GW from an accreting neutron star would provide a unique probe of the 
physics of its interior. 

If the mode is saturated (a > a max ), da/dt = and 

m ^ <^Q + Jil*L p . (3) 



dt r grav 1 - a 2 nSLX Q V 3 / M 

Arras et al. (2002) provide a simple expression for a max in terms of the total energy in the 
r— mode when nonlinear couplings become important, a max = A/(f2|r grav | J). The dimen- 
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sionless constant A (a e in Arras et al. 2002) reflects an inaccuracy in the matching conditions 
between the coupled modes and may range from ~ 10~ 3 to ~ 1. 

A relation determining the thermal evolution of the system closes the system of equa- 
tions. Neutrino emission by the modified URCA process (Shapiro & Teukolsky 1983) is 
taken to dominate the cooling while the core is heated through accretion and dissipation. 
The equilibrium temperature of the core under accretion alone is taken to be T = 10 8 K. 
These assumptions yield the following equation 



dT 

~dt = C, 



V 



a 2 tt 2 MR 2 J n . 8 - 8 , , 
7 x 10 31 (T| - 7f ) ergs s" 1 



(4) 



where C v is the heat capacity of the star taken to be 1.4 x 10 8 T 8 ergs K _1 (Shapiro & 
Teukolsky 1983) and J = 1.635 x 10~ 2 for the density profile considered (Lindblom et al. 
1998). The mass of the star M = 1.4M and its radius R = 12.53 km. 

The calculation begins with T 8 = 1, Cl = 0.1 and a — 0, whenever the system becomes 
unstable to growing r— modes, a is set to 10 -8 . Figure 1 depicts the limit cycle for several 
values of A which determines the amplitude at which the mode saturates (for each value 
of A the peak amplitude of the mode is typically 2 x 10~ 3 A 1 / 2 ). Three important trends 
are apparent. The duration of the GW emission increases dramatically as the saturation 
amplitude decreases. The time for the system to complete the circuit depends only weakly 
on A (3.7-4.8 Myr), so the duty cycle for GW emission also increases. The range of spin 
frequencies over the cycle decreases with A (Andersson et al. 2000). For this model the range 
in spin frequencies over the population of LMXBs also depends on the assumed equilibrium 
temperature of the core T. Stars which accrete less on average will have smaller typical core 
temperatures, and therefore experience a wider range of spin frequencies over the limit cycle. 
Finally, if the r— mode saturates at a smaller amplitude, the temperature of the core during 
the epoch of GW emission decreases. However, for A > 0.1, the temperature of the core 
exceeds 5 x 10 8 K during spin down. 

The relation is well fit by a power law with r on ps 5000A~ a8 yr, in agreement with the 
discussion after Eq. 91 of Arras et al. (2002). The fraction of the time that a source is 
emitting GW increases from ~ 10 -3 to ~ 10 _1 as A decreases from unity to 10~ 3 . To 
probe the possible consequences of a crust, superfluidity or other sources of additional shear 
viscosity, the strength of the shear viscosity is increased by a factor of sixty-four. Here the 
critical spin frequency where the r— mode begins to grow is doubled, the spin-down timescale 
decreases by a factor of 375, and the peak r— mode amplitude increases by a factor of 5.66. 
Decreasing the equilibrium temperature of the core by a factor of two (without changing the 
shear viscosity) has a less dramatic effect of halving the spin-down timescale. Depending on 
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Fig. 1. — The limit cycle in core temperature and spin frequency for an accreting neutron star 
with several values of the dimensionless matching constant, A. The star follows a clockwise 
trajectory in phase space; from the outermost loop inward the values of A are 1, 0.1, 0.01 
and 10~ 3 . The curves are labelled with the duration of the GW emission during each circuit 
of the limit cycle. The bold curve marks the boundary of the instability region. Above this 
curve, the amplification of the r— modes due to gravitational wave emission dominates over 
their attenuation by viscosity. 
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the properties of the neutron star and how the r— mode amplitude saturates, the spin-down 
timescale may vary from 14 to 10 6 years; 34,000 years is the result for a fluid star using 
A = 0.1 (the conservative estimate of Arras et al. 2002). 



2.2. Distribution of Sources - Galactic LMXBs 

Obtaining a better estimate of the number of detectable GW sources requires a estimate 
of the total number LMXBs in the Galaxy, their accretion rates and distances. Since the 
known LMXBs number only about one hundred (Liu et al. 2001), one would expect not to 
observe GW emission from any of them for duty cycles less than one percent. Since LMXBs 
are thought to be the progenitors of millisecond radio pulsars (MSPs), the demographics of 
MSPs can provide an estimate of the total number of Galactic LMXBs (Alpar et al. 1982; 
Bhattacharya & van den Heuvel 1991). Equating the rate of MSP formation to the rate of 
LMXB formation yields an estimate of the number of LMXBs that could be emitting GW 
at any given time, 

AT T ° n AT Ton ( 4- \ — ^MSP r on / k n 
™LMXB,on — iv LMXB — V MSP'-LMXB ) — T on 7 ~ r MSP ~ T77a K* 3 ) 

r C y dc r cycle Afi 10 4 yr 

where r cyc i e is the time for the neutron star to traverse the limit cycle once. This is essentially 
equal to the time for the star to spin up by AQ from the bottom of the limit cycle to the 
top. £lmxb is the time it takes the neutron star to spin up as an accreting LMXB to the final 
spin rate ^msp> an d ?"msp is the rate of formation of MSPs in the Galaxy ~ 3 x 10~ 5 yr" 1 
(Lorimer et al. 1995; Lyne et al. 1998). The ratio of r cyc i c to £lmxb is inferred from Eq. 1. 

Estimating the typical strain amplitude that would be observed from one of these sources 
is also straightforward since the total angular momentum radiated is proportional to AQ. 
This yields 

aR ( 10 1^LGM\ 1/2 1Q _ 25 f_M_ 10VX ^ lOkpc 

8 d \ <n>cH oa ) V1.4M r on J 10km d y) 

where a is a constant which depends on geometry - its mean value is 2.9 (Brady et al. 
1998). To obtain the final approximation, Afi = 0.12 and < VI >= 0.2 which is appropriate 
for A = 0.1. This is thirty times larger than the strain for the brightest source, Sco X-l, 
if the GW emission is continuous (Bildsten 1998) and should be easily detected with the 
initial LIGO (Brady et al. 1998). One could possibly detect such sources even in M31, so 
for r on ~ 10 4 yr, one could expect to find several sources in the Local Group. 

If r on is greater than a few hundred years, the radiation from the surface will reflect 
the heightened core temperature (Tsuruta 1998). The X-ray radiation is powered by the 
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dissipation of the r— modes. The typical core temperature as the star is spinning down is 
6 x 10 8 K, yielding a effective temperature of 3 x 10 6 K and a soft-X-ray luminosity of 
10 35 erg/s (Heyl & Hernquist 2001). Assuming that the source is found at the maximum 
possible distance for an enhanced LIGO with h min ~ 10~ 27 (Brady et al. 1998), 

io~ 27 /ioVV /2 ( M V /2 R 



d - = 7MpC l^-{—) [TamJ ibkm (7) 



this yields an observed X-ray flux of 
F = 2 x 10" 17 erg cnrV 1 
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due to the dissipation of the r— mode energy alone. The typical Galactic source (at a distance 
of 10 kpc), discussed earlier would be 500,000 times brighter with a flux of about 3.2 mCrab. 
Brown & Ushomirsky (2000) discussed the X-ray emission produced if the GW emission is 
steady and derive constraints on either the accretion or the r— modes. If the GW emission 
is transient, the X-ray emission is significantly stronger while GW are being emitted but the 
vast majority of LMXBs would not be in this stage, thereby avoiding these constraints. 



2.3. Distribution of Sources - Extragalactic LMXBs 

For much smaller values of r on , the expected number of active sources in the Galaxy 
vanishes; however, any sources that are active will be visible well beyond the Local Group. 
The total number of sources brighter than /i min is given by 

4 (a R V( m iMl GmY' 2 -i/2»M3P^ , a , 

jYlmxb » = 3* U h~) [ U> 71^—) T - ^T Kmsp (9) 

where IZmsp is the formation rate of MSPs per unit volume averaged over the local region 
of the universe, IZmsp ~ 7"MSp£/-^Miik y Way where C is the local luminosity density, about 
2 x 10 _2 /iiooLMiik y Way Mpc~ 3 (Heyl et al. 1997) where h is the Hubble Constant in units of 
100 km/s/Mpc. Using the values for a max = 1, the case examined by Levin (1999), yields, 

, /10- 27 \ 3 /lyrV 1/2 / M \ 3/ V R V 

The typical distance of the sources detected at this sensitivity with r on ~ 1 yr would be 
nearly 1 Gpc. Even if -/Vlmxb.oii is less than one, one would expect to discover one source 
during a period of r on /iV LM xB,on- 
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3. Discussion 

If the spin of accreting neutron stars is indeed limited by the emission of gravitational 
radiation (Andersson et al. 1999; Bildsten 1998), low- mass X-ray binaries may be an im- 
portant source for LIGO. Although neutron stars may execute a duty cycle (Levin 1999) 
of spin-up and spin-down which reduces the number of active sources at a given time, the 
sources that are active are typically much brighter than in a model where they emit con- 
stantly. Unfortunately, the number of LMXBs that have been discovered actively accreting is 
too small (~ 100) to determine which effect dominates for duty cycles less than ten percent. 

However, if low-mass X-ray binaries are assumed to be the exclusive progenitors of mil- 
lisecond pulsars (see van den Heuvel 1984, for an alternative), the demographics of the 
millisecond pulsars in the Galaxy and beyond provides an estimate of the number of sources. 
Specifically, if the duration of the epochs when the neutron star is emitting gravitational ra- 
diation is greater than 10,000 years, several objects in the Galaxy will be above the detection 
thresholds of LIGO. For r on ~ 10,000 year, these sources would be detectable throughout the 
Local Group. Those objects in the Galaxy could also be detected from their X-ray emission 
which would be powered by gravitational radiation reaction; they would be GW-powered 
neutron stars. 

For r on much less than 10,000 years and much greater than one year, no sources are 
likely to be detectable even with an enhanced LIGO detector. However, if the duration 
of gravitational wave emission per cycle is less than several years, several sources could be 
detected by an enhanced LIGO. In this case, the sources will be located at a typical distance 
of 1 Gpc. 

Depending on the nature of the r— mode instability in the cores of quickly spinning 
neutron stars, LMXBs may provide gravitational-wave beacons throughout the Galaxy and 
the Local Group or at cosmologically significant distances. 

I would like to thank Phil Arras for useful discussions, and I was supported by the Chan- 
dra Postdoctoral Fellowship Award #PF0-10015 issued by the Chandra X-ray Observatory 
Center, which is operated by the Smithsonian Astrophysical Observatory for and on behalf 
of NASA under contract NAS8-39073. 
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